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Abstract 

The Glauber model on a one-dimensional lattice with boundaries (for 
the ferromagnetic- and anti-ferromagnetic case) is considered. The large- 
time behaviour of the one-point function is studied. It is shown that, for 
any positive temperature, the system shows a dynamical phase transition. 
The dynamical phase transition is controlled by the rate of spin flip at the 
boundaries, and is a discontiuous change of the derivative of the relaxation 
time towards the stationary configuration. 
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1 Introduction 



The principles of equilibrium statistical mechanics are well established. But, 
thermal equilibrium is a special case, and little is known about the proper- 
ties of systems not in equilibrium, for example about the relaxation toward the 
stationary state. Some interesting problems in non-equilibrium systems are non- 
equilibrium phase transitions described by phenomenological rate equations, and 
the way the system relaxes to its steady state. As mean-field techniques, gen- 
erally, do not give correct results for low-dimensional systems, people are moti- 
vated to study exactly-solvable stochastic models in low dimensions. Moreover, 
solving one-dimensional systems should in principle be easier. Exact results for 
some models on a one-dimensional lattice have been obtained, for example in 
[1-14]. Different methods have been used to study these models, including ana- 
lytical and asymptotic methods, mean field methods, and large-scale numerical 
methods. 

The Glauber dynamics was originally proposed to study the relaxation of 
the Ising model near equilibrium states. It was also shown that, there is a re- 
lation between the kinetic Ising model at zero temperature and the diffusion 
annihilation model in one dimension. There is an equivalence between domain 
walls in the Ising model and particles in the diffusion annihilation model. Ki- 
netic generalizations of the Ising model, for example the Glauber model or the 
Kawasaki model, are phenomenological models and have been studied exten- 
sively [15-20 ]. Combination of the Glauber and the Kawasaki dynamics has 
been also considered [21-23]. 

In , an asymmetric generalization of the zero-temperature Glauber model 
on a lattice with boundaries was introduced. It was shown there that, in the 
thermodynamic limit, when the lattice becomes infinite, the system shows two 
kinds of phase transitions. One of these is a static phase transition, the other 
a dynamic one. The static phase transition is controlled by the reaction rates, 
and is a discontinuous change of the behavior of the derivative of the station- 
ary magnetization at the end points, with respect to the reaction rates. The 
dynamic phase transition is controlled by the spin flip rates of the particles at 
the end points, and is a discontinuous change of the relaxation time towards the 
stationary configuration. Other generalizations of the Glauber model consist 
of, for example, alternating-isotopic chains and alternating-bound chains (see 
[ p5| , for example). People have also considered phase transitions induced by 
boundary conditions (see [26-28], for example). 

The scheme of the paper is as follows. In section 2, the model is introduced, 
the rates are determined using the detailed balance, and the steady state con- 
figuration of the magnetization is obtained. In section 3, the dynamical phase 
transition of the system is investigated, and it is shown that it does show a 
dynamical phase transition, provided the temperature of the system is not zero. 
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2 Kinetic Ising model on a one-dimensional lat- 
tice with boundaries 



The model being addressed, is the Glauber model on a one-dimensional lattice 
with boundaries. In the Glauber model, the interaction is between three neigh- 
boring sites. Spin flip brings the system into equilibrium with a heat bath at 
temperature T. A spin is flipped with a rate fi = 1 — tanh(2/3J) whenever the 
spin of both of its neighboring sites are the same as it and is flipped with a 
rate A = 1 -I- tanh(2/3 J) whenever the spin of both of its neighboring sites are in 
the opposite direction. (Here f3 := 1/{KT).) At domain boundaries, spins are 
flipped with unit rate. So the interactions can be written as, 

TTT^Tit andiii^iti with rate ^ 
Tit^TTT and iti^iiiwithrate A 
T T i ^ T i i and i i T ^ i T T with rate 1 

Consider a lattice with L sites and the Glauber dynamics as the interaction. 
The spin of the first site may flip with the following rates 

T i iwith rate gi 

T i Twith rate g2 

i T Twith rate 53 

i T jwith rate gi, 

and the spin of the last site may flip with the following rates 

i i iwith rate hi 

T T iwith rate /i2 
T i~* T Twith rate 

i i^ i Twith rate ^.4. 

It is known that the time evolution equation for the one-point functions in the 
bulk are expressed in terms of only the one-point functions To make this 
true for the boundaries as well, the following relations should hold. 

91 + 54 =52 + .93 

hi + hi ^h2 + hs. (1) 

One may give a physical meaning to the parameters gi and hi , by demanding 
the detailed balance to hold. Consider the energy £ of the system to be 

L-l 

f = -fSiSL-f J ^SiS,;+l), (2) 
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then, the detailed balance demands 



7?-(siS2 ^ s[s2) exp{/3(BiSi + JsiS2 H )} 

=7^(s'lS2 ^ S1S2) exp{p{Bis[ + Js[s2 + ■■■)}, (3) 

where TZ{siS2 — > S1S2) is the rate of the spin flip of the first site from si to s'l- 
Equation (||) shows that 

7^(slS2 s[s2) = f{s2) exp{-/3si(Bi + JS2)}. (4) 

The exponential term in the above equation is at most linear in terms of si. So, 

TZ{siS2 s[s2) = /(S2)[l - si tanh/3(Si + JS2)]. (5) 

Then, 

5i = /(-l)[l-tanh/3(Bi- J)] 
52-/(l)[l-tanh/3(Bi + J)] 

53 = /(l)[l + tanh/3(Bi + J)] 

54 = + tanh/3(Bi - J)]. (6) 

The condition of exact solvability (|l|) (the closure of time evolution equation of 
one-point functions) leads to 

/(I) = /(-!)■ (7) 

This means that the inertia of the first spin against spin-flip does not depend on 
the second spin. A similar expression can be written for the rate of the spin-flip 
of the last site. 

For the infinite lattice, the Glauber model has a particle reaction-diffusion 
interpretation. If the spins of the neighboring sites are different (at a domain 
wall), one may consider the link between that sites as a particle. When the 
spins of the neighboring sites are the same (no domain wall) , one may consider 
the link between the sites as a vacancy. Then the Glauber model turns into a 
reaction-diffusion model: 

• • ^ o owith rate f + tanh(2/3 J) 
o o ^ • •with rate f — tanh(2/3 J) 

• 0^0 •with rate 1, 

where a particle (a vacancy) is denoted by • (o). For the Glauber model with 
boundaries, to have a consistent particle model, one has to impose 

gi =53 52 = 54 

hi = hs /i2 = Ih. (8) 
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Then, the injection and extraction of particle at the first site are 

• owith rate gi = 

o •with rate (72 =94, (9) 

and the injection and extraction of particles at the last site are 

• — > owith rate hi — /13 

o — !• •with rate /12 = h^, (10) 

Now, consider the general case where only the conditions (|^), which guarantee 
the closure of the time evolution, are satisfied. We have 

(sfe) = - 2{sk) + ((sfe+i) + (sk-i)) tanh(2/3J) l<k<L 

(si) =-{92+ 93){si) + (51 - 52) (S2) + (33 - gi) 

(sl) =-ih2 + h3){sL) + (hi - h2){sL-i) + {hz - h^). (11) 

The steady-state solution to (^Tj) is 

{sk)=Diz'l + D2Z^-'^-\ (12) 

where 

zi = 22^1 = tanh(/3J). (13) 
It can be shown that in the thermodynamic limit [L 00), 

91 -93 



Di = 
D2 



{91 - g2)zl - {92 + 93)21 

hi - /13 



{hi - h2)zf - (/i2 + h3)zi ' 

Di, and D2 are continuous functions of the rates. So the behavior of (s^) near 
the ends of the lattice varies continuously with rates, and there is no phase 
transition. 

3 The dynamical phase transition of the system 

The average magnetization per site m{t) is 

^W = |E(^;.W)- (14) 

fc=i 

In the thermodynamic limit, the boundary terms are negligible, and 

-^m(i) = 2[tanh(2/3J) - l]m{t). (15) 
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Then, the same as Glauber model on an infinite lattice, the average magneti- 
zation do not show any phase transition. But, as it will be shown, the system 
does exhibit dynamical phase transition. 

The homogeneous part of ( pi] ) can be written as 

{h) = hUsi). (16) 

The eigenvalues and eigenvectors of the operator h satisfy 

E Xk^ -2xk + tanh(2/3J)(a::fc+i + Xk-i), k^l,L 
E xi ^ -{g2 + g3)xi + {gi - g2)x2, 

E XL ^ ~{h2 + h3)xL + {hi - h2)xL-i, (17) 

where the eigenvalue and the eigenvector have been denoted by E and x, re- 
spectively. The solution to these is 

Xk=az'l+bz^, (18) 

where 

-{E + .92 + g3.)iazi + bz2) + {gi - g2){azl + bzl) = 
-{E + h2 + h3){az^ + bz^) + {hi - /i2)(azf + bz^-^) = 0, (19) 

and Zj's satisfy 

£; = -2 + tanh(2/3J)(^ + ^"i). (20) 

So, ziZ2 — 1. Two cases may occur, 

a) zi and Z2 are phases. 

b) zi and Z2 are both real, not equal to ±1. Then the modulus of one of them 
is less than one, that of the other is greater than one. 

Using ( po|) and ziZ2 = 1, one can eliminate E, and arrives at 

7?^^[2 - 52 - .93 + z{gi - .92 - tanh(2/3J)) - z^^ tanh(2/3J)] 
X [2 - /i2 - + z{hi - /i2 - tanh(2/3J)) - z^^ tanh(2/3J)] 
-z^-^[2 - 32 - 33 + z^^igi - .92 - tanh(2/3J)) - ztanh(2/3J)] 
X [2 - /i2 - /i3 + z-^{hi -h2- tanh(2/3J)) - z tanh(2/3J)] = 0. (21) 

Obviously, Zj = ±1 satisfies (pl|). But these solutions lead to 

Xk = z^{a + bk). (22) 

And this form for x/c, generally does not satisfy the boundary conditions at 
fc = 1,L. Equation (|l|) can be written in the form 

G{z):^F{z)~F{z-^) = 0, (23) 
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where 



F{z) z^-^{2 - g2 -93 + z[gi -92- tanh(2/3J)] - z"^ tanh(2/3J)} 

x{2 - /i2 - /i3 + z[hi -h2- tanh(2/3J)] - z"4anh(2/3J)}. (24) 

For a phase solution to (21 ), 2; = e*'', we have 

£; = -2 + 2 tanh(2/3 J) cos I?. (25) 

In the thermodynamic limit (L —>■ go), in any neighborhood of z = 1 there 
exist a phase solution to (|2l] ) . The supermum of the eigenvalues determines the 
relaxation time toward the stationary average-density profile. So, if all of the 
solutions are phase, 

T = [-2 + 2 tanh(2/3J)]-^ (26) 

But, if there exist solutions which are not phases, they should be real. Consider 
z > 1. Then for L 00, ( pT] ) becomes 

[2 - 32 - 33+^"^ (.91 - .92 - tanh(2/3J)) - z tanh(2/3 J)] 
2 - /i2 - /i3+z"i(/ii - /i2 - tanh(2/3J)) - z tanh(2/3 J)] = 0. (27) 



X 



Changing the rates, one may arrive at a situation where the above equation has 
a real solution greater than one. The transition occurs at the point that this 
equation has a solution equal to one. When the system has passed this point, 
the relaxation time becomes 

T = [-2 + 2(A + A-i) tanh(2/3J)]-\ (28) 



where A is that solution to (]27|), which is greater than one. (here we have 
assumed J > 0, ferromagnetism. If J < 0, anti-ferromagnetism, A is that 
solution to ( |27| ) which is less than —1.) Putting z = 1 in (p7|), at least one of 
the following equations should hold 

2[1 - tanh(2/3J)] - 52 - .94 = 

2[1 - tanh(2/3J)] - /i2 - /i4 = (29) 

If the temperature is zero, (^9|) for example gives 52 +.94 — 0. Remembering that 
these parameters are rates, one arrives at 92 = .94 = 0. So, at zero temperature, 
the solution cannot pass z = 1. But at any other temperature, 1 — tanh(2/3J) is 
positive, and changing the parameters, 92 + 54 can be made more than or less 
than 1 - tanh(2/3J). 

If one uses the expressions (^) and ^ for g^'s, then ( |29| ) becomes 

-2tanh(2/3J)+/[tanh/3(J-Bi)+tanh/3(J + Bi)] ==0. (30) 

Putting / = 1, means that the inertia of the first spin against the spin flip is 
the same as those of the bulk spins. In this case, however, (|30|) has no solution. 
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That is, there is no phase transition. In fact, ( pO[ ) has no solution for / < 1. 
For / > 1, however, it may have a solution. 

It is seen that the parameters 32 and (74 (or /12 and /14) are control parameters 
of the dynamical phase transition. The parameters gi and 53 (or hi and /13) do 
not have any contribution in the dynamical phase transition. The rates gi and 
g3 are the rates of the disappearance of the domain walls. But we note that the 
eigenvector corresponding to z = 1 is a configuration where all the spins are the 
same {sk ~ z*^ = 1.) It is this configuration which corresponds to the largest 
value of E, which determines the relaxation time, and in this configuration, there 
is no domain wall. The disappearance rate of this configuration determines 
the relaxation time towards the steady state, and gi and 173 (or hi and /13) 
are irrelevant to this rate. In the particle-vacancy picture, this means that 
the rate of change of vacancy to particle is important, since the configuration 
corresponding to the maximum value of E is the empty lattice. 

This arguments are true for J > 0, the ferromagnetic case. If J < 0, then 
the relaxation time is determined by the value of E at the smallest possible 
value of z (which is less than —1), and the transition occurs as z = — 1 becomes 
a solution to (|2^). It is not difficult to see that in this case gi + g^ (or hi + /13) 
determine the phase transition. The reasoning is the same as above, except that 
here the configuration determining the relaxation time is that corresponding to 
z = — 1, which means that the spins are alternating. So, in this configuration 
there are no ]] or [[ configurations and g2 and 54 (or /12 and /14) are irrelevant. 

Acknowledgement The authors would like to thank Institute for Studies in 
Theoretical Physics and Mathematics for partial support. 
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